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SUMMARY 


hiKTgy  solutions  are  excellent  analysis  procedures  lor  predicting  residual  strains  or  de- 
formations in  structural  elements  when  transient  behavior  is  of  little  interest.  In  this  report, 
we  demonstrate  how  energy  procedures  can  be  used  in  rigid-plastic  structural  solutions  when 
memhersarc  loaded  either  impulsively  or  quasi-stalically  by  blast  waves. 

In  the  impulsive  loading  realm,  the  kinetic  energy  imparted  to  a structural  member  is 
equated  to  the  plastic  strain  energy,  whereas  in  the  quasi-statie  loading  lealm.  the  work  per- 
formed in  deforming  a structural  member  is  equated  to  the  plastic  strain  energy.  An  assumed 
first  mode  structural  deformation  pattern  works  well  when  calculating  plastic  strain  energy 
in  either  of  liie  loading  realms.  Experimental  test  data  on  defonned  simply-supported  and 
anliiever  beams,  clamped  eircidar  plates,  and  clamped  rectangular  plates  demonstrate  the 
validity  of  these  solutions. 

The  test  data  on  different  types  of  structural  elements  are  important,  as  the  beam  data 
involve  only  bending  behavior,  the  circular  plates  have  both  bending  and  extensional  action, 
and  the  rectangular  plates  introduce  .shearing  behavior  into  the  strain  energy  calculations. 
Because  all  solutiojis  are  clo.sed-form  ones,  design  fomnilae  result  which  can  be  used  to  eval- 
uate plastic  dcfornuition  in  blast  loaded  structural  members. 

This  report  is  a reprint  of  a paper  presented  at  the  Iblh  Explosive  Safety  Seminar. 
Hollywood  Beach.  Florida,  September  |d74. 
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ENER6Y  SOLUTIONS  PGR  PREDICTING  DEFORMATIONS 
IN  BLAST  LOAOEO  STRUCTURES 


I.  INTRODUCTION 

Energy  solutions  arc  excellent  analysis  procedures  for  predicting  residual  strains  or 
deformations  in  structural  components  whenever  transient  (time-dependent)  behavior  is  of 
little  interest.  Although  schools  readily  teach  energy  procedures  for  obtaining  elastic  solu- 
tions, few  investigators  apply  the  approach  to  dynamic  plasticity  rroblcms. 

Dynamic  rigid-plastic  energy  solutions  began  in  the  early  1950’s  when  Lee  and 
SymondsI  * I used  the  static  plastic-hinge  concept,  considered  beam  inertia,  and  propagated 
a traveling  hinge  to  analytically  obtain  '.he  upper  bound  for  permanent  deformation  in  a 
beam  under  a transverse  load.  Their  Brown  University  associates  and  graduates  such  as  S.  R. 
Bodner.  W.  Prager.  N.  Jones.  J.  B.  Martin.  R.  M.  Haythornthwaite  and  others  then  added 
refinements,  illustrations  of  which  arc  given  in  Refs.  2-6.  It  was  J.  E.  Greenspon  in  the  I960's 
who  pointed  out  that  one  could  obtain  solutions  without  going  through  the  details  of  propa- 
gating a plastic  hinge  along  structural  members  I Refs.  7-111.  Greenspon  noted  that  the  resid- 
ual strain  energy  stored  in  a plastically  deformed  member  could  be  calculated  by  assuming  a 
final  deformed  shape.  This  strain  energy  was  then  equated  to  the  energy  flux  in  an  explosive 
blast  wave. I * ’ I We  disagree  with  this  last  step,  which  made  deformations  independent  of 
structural  orientation  relative  to  the  enveloping  blast  wave,  thus  ignoring  an  important  effect 
observed  in  many  e.xperimonts.  In  addition.  Greenspon’s  procedure  forces  pressures  and 
impulses  in  the  blast  wave  to  obey  the  relationship  PI  = constant,  whose  asymptotes  for  both 
pressure  and  impulse  are  Z’  = 0 and  / = 0.  Tlie  response  of  real  targets  is  related  to  non-zero 
P and  / limits,  so  this  conclusion  is  also  unacceptable. 

We  would  emphasize  that  Greenspon  was  correct  in  his  strain  energy  estimation  procedures. 
Estimates  of  structural  deformation  would  have  been  correct  had  he  equated  strain  energy  to 
the  kinetic  energy  imparted  to  the  structure  for  short  duration  impulsive  loads.  When  durations 
are  long  relative  to  the  structural  response  time,  the  strain  energx'  is  equated  to  the  work  per- 
formed when  the  peak  load  moves  through  the  distance  that  the  structure  deforms.  Hence,  two 
separate  procedures  are  required,  one  to  obtain  the  solution  for  the  impulsive  loading  realm, 
and  the  other  to  obtain  the  solution  for  the  quasi-static  loading  realm.  We  will  proceed  to  illus- 
trate these  procedures  by  computing  results  and  comparing  the  test  data.  Our  first  illustration 
is  a rheological  model  whose  exact  solution  can  be  obtained  and  compared  to  the  answers  given 
by  energy  procedures. 

II.  SINGLE  DEGREE-OF-FREEDOM  SYSTEMS 

Consider  first  a single-degree-of-frccdom.  rigid-plastic  system  as  in  Figure  la.  The  motion 
of  the  mass  m is  resisted  by  a Coulomb  friction  element  / when  the  blast  load  /)(/)  is  applied  to 
the  structure.  Wc  will  appro.ximatc  the  blast  loading  with  an  exponential  decay  as  in  Figure  lb 


r 
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HCURb  I.  HIGID-l'LASTIC.SINGLb-DECRbE-OI -FREbDOM  DAMAGE  MODEL 


(whcie  Z’  is  the  maximum  applied  force  and  T is  the  time  constant  associated  with  the  duration 
oi  loading).  M PI)  ^ 1 .0.  wc  have  the  trivial  case  where  the  residual  deformation  X equals  0 
because  the  mass  never  moves.  If  P!f>  1.0,  we  can  write  the  differential  equation  of  motion: 

ft’ ") 


By  direct  integration,  we  obtain  for  the  ca.se  of  zero  initial  velocity,  the  velocity  relationship: 

(2) 


d/  L \p/\r)_ 


Integrating  again,  we  obtain  for  the  case  of  zero  initial  displacement,  the  displacement  equation 

PP 


X = ■ 


in 


(3) 


Motion  continues  until  the  velocity.  Eq.  (2),  equals  zero  or  until: 


(4) 


We  cannot  explicitly  solve  for//?' in  Eq.  (4).  as  it  is  a transcendental  equation:  therefore, 
we  assume  values  of  PI/,  solve  for  tiT.  and  substitute  into  the  displacement  equation  (Eq.  (3)1 
to  obtain  the  maximum  deformation  X.  Table  1 gives  the  results  of  such  a calculation. 

The  maximum  deformation  X in  the  third  column  of  Table  I has  been  made  nondimen- 
sional  by  dividing  the  left  and  right  sides  of  Eq.  (3)  by  (Pp  )lin.  A solution  can  be  presented 
for  the  maximum  displacement  by  plotting  (Xni)l(PT- ) versus  Pif.  We  have  elected  to  divide 
PI/ by  (Xm)l\PT^)  to  form  a new  fourth  column  in  Table  I and  to  plot  this  new  column 
{PP^  l(XmJ)  versus  P!/  The  rexson  for  this  manipulation  is  that  the  product  PT  equals  the 
applied  total  impulse  /,  and  in  this  manner  we  create  a scaled  load-impulse  or/'-/  diagram. 

The  solid  line  in  Figure  2 is  this  scaled  P-l  diagram  for  a simple  rigid-plastic  structure. 
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Ohscrs'c  in  Tabic  1 and  Figure  2 lhat  wlicncvvr  /•  liXinf)  is  greater  than  about  60.  the 
duration  of  loading  Tis  larger  than  the  response  tinie  / and  /Y/ equals  1.0.  Similarly,  when- 
ever/V/ is  greater  than  about  20.  durations  of  loadwig  7'are  smaller  than  the  response  limes  / 


IIMMI • ^ 


of  the  respondini  structure,  and  /*  liXmf)  equals  2.0.  The  enerfy  solutions  which  we  will 

apply  estimate  both  of  these  ^ymptotcs. 


The  strain  energy  U stored  ir.  plastic  deformation  is  given  by: 

U^fX 

The  kinetic  energy  KE  imparted  to  the  mass  equals; 


(5) 


t6) 


The  work  h'  done  by  the  maximum  force  P acting  through  the  distance  X is; 

h’  = PX 


(7) 


Equating  the  strain  energy  U.  Eq.  (5),  to  the  kinetic  energy  KE,  Eq.  (6),  yields  the  asymptote 
for  the  impulsive  loading  realm  or: 


Xmf 


= 2.0 


(impulsive  loading  realm) 


(8) 


The  other  asymptote  is  obtained  by  equating  U.  Eq.  (5),  to  the  work  lE,  Eq.  (7). 

(quasi-static  loading  realm) 


J 


(9) 


Had  we  wished  to  calculate  maximum  clastic  deformations  rather  than  plastic  ones,  the 
same  procedures  would  apply.  Replacing  the  Coulomb  friction  element  with  a linear  elastic 
spring  in  Figure  la  would  have  yielded  the  analytical  solid  curved  line  shown  in  Figur-.-  1. 

Tliis  solution  also  has  asymptotes  for  the  impulsive  and  quasi-static  loading  realms  that  can  be 
obtained  U'^ing  energy  procedures.  The  strain  energy  in  an  elastic  system  would  be  given  by 
Eq.  ( 1 0)  rather  than  Eq.  (5). 
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U = ^kX^ 


(10) 


In  an  elastic  system  the  kinetic  energy  KE  and  work  h'arc  still  given  by  Eqs.  (6)  and  (7). 
respectively.  For  an  elastic  system,  equating  Eq.  (10)  to  Eq.  (6)  yields  the  asymptote  for  the 
impulsive  loading  realm. 


or 


(impulsive  loading  realm) 


(12) 


A. 

p/k 


'~T 

ni 


Similarly,  equaling  tq.  ( 10)  to  Hq.  (7)  yield.s  the  quasi*statie  loading  realm  asymptoie; 


(quasi-slalie  loading  realm) 


(i.^) 


Both  the  impulsive  loading  realm  asymptote  and  quasi-stalie  loading  malm  asymptote  are 
shown  in  Figure  3.  Figure  3 illustrates  that,  if  these  asymptotes  are  known  lor  an  elaslie  as 
well  as  plastic  system,  the  deformations  can  be  predicted  using  energy  pr<Kediires. 

As  has  been  illustrated  in  these  simple  models,  the  principles  are  as  follows; 

1 1 ) To  estimate  the  impulsive  loading  realm  strueiural  defomialion  asymptote,  estimate 
the  strain  energy  in  a deformed  .structure  end  equate  this  strain  energy  to  the  ki.netie 
energy  imparted  to  the  structure. 

(2)  To  obtain  the  riuasi-static  structural  uetormatiim  asymptote,  eijua'e  the  strain  energy 
to  the  work  performed  by  the  peak  force  deneeling  with  the  siruciure. 

We  are  now  prepared  to  illustrate  these  prinei|>Ies  as  they  are  applied  to  beam,  plates,  and  simi- 
lar more  complex  structural  components  which  are  loaded  by  blast  waves  or  other  transient 
pulse..  The  Hr.st  structural  components  that  we  will  study  are  cantilever  beams. 

III.  BLAST  LOADED  CANTILEVER  BEAMS 

Consider  a damped  cantilever  beam  of  rectangular  cross  section  as  dejucied  in  l igure  4. 
We  will  a.ssumc  that  the  defomied  shape  for  this  structural  component  is  given  by 


Notice  that  the  assumed  defomied  shape  has  no  detlection  and  no  slope  at  v = 0 and  no 
moment  at  \ = L.  The  maximum  deformation  and  maximum  slope  occur  at  v = /..  a:id  the 
ma.ximum  moment  occurs  at  .v  = 0.  These  are  the  correct  boundary  conditions.  .*\  good  solu- 
tion does  depend  upon  selecting  an  appropriate  defomied  shape. 

Because’  no  membrane  action  is  developed,  the  strain  energy  is  d'ssipated  in  bending.  In 
an  elastic  member,  the  bending  strain  energy  is  given  by: 


I I 


riGURE  3.  SfAXlMUM  RKSPONSK  TO  i ORCt  PU.'.Si: 


I KU'RK  4.  sail  MATIC  DRAWINc;  OI  fl  AMPI  D 
CANIII  I VI  k hi  am 


.1: 

whiTv  (he  moment  M is  equal  to  /;7(d^  *e/tlA'’ t.  DilTerentiating  Kq.  ( 14)  and  Mihstiliitin^ 
into  iiq.  (15)  then  yields: 


rr'’/:7uv,  / , /7Ta  \ , 

““  (57)*'' 


Or,  upon  integrating: 


IT  IHw;, 

li  = ^ — ( bending  in  elastic  cantilever  beam ) ) 

ML 

The  kiiK  tic  energy  imparled  to  the  cantilever  beam  <'!’  mass  density  p.  width  />.  and 
thickness  /;  by  a iinifomi  specific  impulse  of  intensity  / is  given  by: 

A7.  = 23  ~ tt/  Tii  ( 


r I / ih  V \‘ 

A A'  = f -(phh  d \ )( ) 

“ : VpA/zdv/ 


Integration  of  the  preceding  quantity  then  yields: 


r 


(20) 

2ph 

The  equation  deOning  the  impulsive  loading  realm  asymptote  for  the  elastic  maximum  ti;/ 
deformation  is  then  obtained  by  equating  Eq.  (20)  to  Eq.  < 1 7).  The  width  b drops  out  of  the 
solution  when  M^/12  is  substituted  for  the  second  moment  of  area  / in  Eq.  (17). 


^ / / \ 
L \/t)  \h\/Ep) 


The  quasi-static  loading  realm  asymptote  is  pn.*dicted  by  estimating  the  work  associ- 
ated with  the  peak  drag  load  Q.  The  work  is; 

122) 


’ = (-'o  Qbn'u  / ('  “ 


Integration  of  tlic  preceding  quantity  then  yietds: 


Equating  Eq.  (24)  to  Eq.  ( 17)  and  substituting  for  the  second  moment  of  area  then  yields  the 
quasi-static  asymptote  for  the  elastic  ma.xiimim  tip  deformation. 


7 = ^(^)(^)’ 


Maximum  elastic  strains  at  the  njot  of  the  cantilever  beam  can  also  be  calculated  from 
the  ma.ximum  tip  deformations,  either  Eq.  (2 1 lor  Eq.  (2.S).  dependent  upon  the  loading  realm. 
Substituting  the  second  derivative  of  Eq.  ( 14)  into  the  moment  curvature  relationship,  the 
moment  into  a - Mc/l.  the  stress  u into  c = «//:'. and  A/rVl2  for  the  second  moment  ■•f  area  /. 
yields  the  strain  eqi;ation  in  terms  of  the  maximum  tip  defonnation.  This  cqua'ion  is  given  by: 

TT^/m’n  ir.v 


The  strain  e will  be  a maximum  at  the  root  of  the  cantilever  beam  where  cos  Tr.\/2L  equals  1 .0. 
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Substituting  Eq.  (2 1 ) for  Hig  into  Eq.  (26)  then  yields  the  maximum  strain  equation  for  the 
impulsive  loading  realm. 

e = 3.45 — (impulsive  loading  realm)  (27) 

hi/tp 

I * 

Equation  (26)  for  the  impulsive  loading  realm  can  be  compared  to  experimental  test 
data  to  demonstrate  the  validity  of  this  analy.Ms  procedure.  In  1958.  Baker,  el  al..l ' - I 
1 detonated  H.E.  explosive  charges  in  the  vicinity  of  606I-T6  aluminum  cantilever  beams. 

Plotted  in  Figure  5 is  the  maximum  bending  strain  as  a function  of  i/Uy/^)  for  beams  with 
a length  of  12  in.  and  thickness  of  0.05 1 inch.  Some  uncertainty  exists  in  computing  the 
impulse  imnarted  to  the  beams  because  of  an  air  blast  wave  diffracting  around  the  beams; 
hence,  the  test  data  are  plotted  as  bars.  As  can  be  seen  in  Figure  5.  Eq.  (27)  for  the  maximum 
elastic  strain  at  the  root  of  the  cantilever  beams  accurately  predicts  experimentally  uUserved 
results. 


'I  he  strain  at  the  root  of  cantilever  beams  in  the  quasi-static  loading  realm  can  be  obtained 
by  substituting  Eq.  (25)  into  Eq.  (26)and  setting  the  cos  I(7rC)/(2Z.)l  equal  to  1.0.  Tliis  result 
yields: 


(quasi-static  realm) 


(28) 


J.  D.  Day  in  unreported  tests  ran  five  blast  loading  experiments  on  clamped  cantilever  steel 
beams  with  strain  gauges  at  the  root  of  the  beam.  The  steel  beams  were  6 X 3/4  X 3/4  in.  and 
were  exposed  to  side-on  overpressures  Pj  as  summarized  in  Table  2.  The  drag  pressure  Q can 
be  calculated  from  P,.  and  assuming  that  Co  equaled  1.75,  we  obtain  calculated  strains  as 
shown  in  the  last  column  of  Table  2.  These  calculated  strains  compare  very  favorably  with 
experimental  observed  results. 

1 

Plastic  response  of  cantilever  beams  can  be  estimated  as  easily  as  clastic  response.  If  we 
assume  a rigid-plastic  model,  the  bending  strain  energy  now  becomes  the  integration  of  the 
moment-curvature  relationships  over  the  length  of  the  beams.  In  other  words: 


/VI  IV 

My  — d.X 

0 


d*»v 


(29) 


Differentiating  Eq.  (14).  the  assumed  deformed  shape  and  substituting  it  into  Eq.  (29)  gives: 


i'  Tt.x 


U 


2 f ft  ' 


(30) 


MAXIMUM  BEMDING  STRAIN  x 10 
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TA«Ui  2.  ELASTIC  STRAINS  IN  CANTILtVKRS. 
QUASI-STATIC  LOADING  REALM 


p,  (psl) 

Q(pit) 

Experimental 
£ X IO+‘ 

Calculated 
e X t0^‘ 

n 

8. HO 

n 

23.2 

23.6 

H.6K 

WSSm 

22.9 

23.4 

10.48 

32..S 

■■ 

16.19 

80.4 

73.5 

iHli 

20.78 

8.7.t 

141.0 

117.5 

But  in  a rectangular  member  the  yield  moment  My  = af>h^l4:  hence,  upon  substituting  for 
/»/,.  and  completing  the  integration,  we  obtain: 


nOyhh^Wf, 

I6i 


(31) 


Kquating  l{q.  (3 1 ) fo  the  kinetic  energy,  tq.  ( 20).  then  yields  the  asymptote  tor  perma- 
nent residual  tip  detlection  in  the  impulsive  loading  malm. 


I'V.  8//.W  / Y 

L ir\///\//V^/ 


(32) 


Reference  I 2 also  presents  tost  data  on  permanent  tip  deneetions  for  the  same  606  l-K) 
aluminum  cantilever  beams  as  were  used  in  the  elastic  strain.  Figure  .S.  comparison.  Plotted 
in  Figure  6 are  these  sealed  tip  deflections  as  a function  of/7(£>//^).  Kqiiation  (32)  was 
placed  in  the  same  tbrmat  as  the  data  in  F'igure  6 bv  multiplying  and  dividing  F'q.  (32)  by  the 
elastic  modulus  E.  As  can  be  seen  in  Figure  o,  the  agreement  is  relatively  good.  The  disagree- 
ment that  does  arise  at  small  values  of  scaled  impulse  //(Z.v'nTT)  is  caused  by  our  use  ol  a tigid- 
plastic  rather  than  elastic-plastic  analysis 

The  writers  have  no  data  for  a cantilever  bcani  permanent  tip  deflection  comparison  in 
the  quasi-static  loading  realm.  In  the  quasi-static  loading  realm,  the  displacement  .V  bcLomes 
indetemiinate  when  V,  F'q.  (3 1 ).  is  equated  to  ll'A'.  Fa|.  ( 24).  The  .solution  in  the  quasi-static 
loading  realm  is  uiven  by: 


P //f\  = 


(33) 


Tile  displacement  ivo  in  F'q.  (33)  is  indeterminate  just  as  X was  indeterminate  in  the  rigid- 
plastic  rheological  model.  Fq.  (9).  for  the  quasi-static  loading  realm.  This  conclusion  is  cor- 
rect for  perfectly  plastic  nonhardening  systems  in  the  quasi-static  loading  realm.  A hardening 
stress-strain  law  dot-s  result  in  finite  displacements. 
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IV.  BENOING  IN  SIMPLY  SUPPORTED  AND  CLAMPED  BEAMS 

Our  next  illustration  will  be  plastic  bending  in  a simply-supported  beam  being  loaded 
with  n uniform  load.  Figure  7 shows  the  deformed  shape  of  a bent  simply-supported  beam. 


L — 

I I i i i i t t 



HCURt  7.  OtFORMtU  SlUPK  OF  SIMPLY  SUl*K)RTtD  BFA.M 
To  calculate  strain  energy  in  this  member,  we  must  assume  a deformed  shape.  Selecting; 


as  an  appropriate  deformed  shape  gives  the  appropriate  boundary  conditions.  At  v = 0. 
center  of  the  beam,  u’  = Uq  and  the  slope  (dw/dv)  = 0,  while  at  the  ends  of  the  beam,  v = 
till,  the  deflections  u-  = 0 and  the  slopes  (dtv/dv)  = ~8.vvi>o/Z,^ . a maximum  value.  The 
strain  energy  equals  the  plastic  yield  momc:it  My  for  the  beam  cross-section  times  the  change 
in  angle  of  rotation  integrated  over  the  entire  beam.  Because  the  beam  is  symmetric  and  the 
change  in  angle  of  rotation  With  respect  to.v  approximately  equals  — d^u  /dv-,  the  strain 
energy  equals; 


i.h 

V = ~2f  My  — dv 
0 


dv 


.2 


(35) 


Differentiating  Eq.  (34).  substituting  it  into  Eq.  (35)  and  integrating  then  yields: 

y _ 8T/|.  tt’o 


The  kinetic  energy  KE  is  obtained  by  summing  the  impulse  squared  divided  by  two  times 
the  incremental  mass  (sec  Eq.  (6)|  over  the  entire  beam.  If  b is  the  width  of  the  loaded  mem- 
ber. p the  density.  A the  cross-sectional  area,  and  / the  specific  impulse,  this  summation  yields 
the  following  integration. 


/ 2pA{d.\) 


(37) 
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or 


I 


I 

t 


IpA 


(38) 


Equating  U.  Eq.  (36).  to  KH,  Eq.  (38),  yields  the  asymptote  Tor  the  impulsive  loading 
realm. 


i^h^L  /‘v'o\ 

— — : = 16(  T 1 (impulsive  realm  s.s.  beam) 

piMyA  \L/ 


(39) 


The  work  W is  obtained  by  integrating  over  the  length  of  the  beam,  the  forces  times  the 
distances  through  which  they  move.  This  operation  is  performed  by  integrating  dv  times 
the  assumed  deformed  shapes.  Eq.  (34). 


r./2  / 4v*\ 

It'  = 2 / phwo  {} 


(40) 


or 


ll' = jpbf,H-o  (41) 

Equating  It'.  Eq.  (41),  to  U,  Eq.  (36).  yields  the  quasi-static  asymptote, 

— — =12  (qiiasi-static  realm  s.s.  beam)  i42) 

My 

So  far  these  calculations  have  assumed  that  the  beam  is  simply-supported  and  free  to 
rotate  at  the  ends.  If  the  beam  is  clamped  at  the  ends,  no  rotation  occurs,  but  can  move 
inwards  so  that  no  membrane  action  is  developed:  we  can  use  many  of  the  results  which  have 
already  been  developed.  To  do  this,  assume  that  a clamped  beam  is  really  two  simply- 
supported  beams  that  have  been  split  and  joined  end  to  end  as  in  Figure  8. 

This  new  contlguration  implies  that: 


u = 


I (>■>/>•  Jo 
V 


and  that 


A/.-  = :ay.\.,.  = 


pA 


(43) 


(44) 


20 


i 


i 


ini'll '^l•'!l•lllIl■l«^^w•'"Pl*W'l 


I ICl'RI  8 Df  l-ORM^:i)  shark  OI  CKAMM-  D BK  AMS 


but  ^ = L '2  and  ,»o  ~ Mo/->  so  substituting  for  (!  and  ,i  o *H)ii  equaling  U to  KE  yields: 


(impulsive  realm  clampcii  beam) 


Because  Kq.  (4S)  for  clamped  beams  is  twice  Eq.  (39)  for  simply-supported  beams,  we 
can  write 


—JllL  = ]/,  ^ 

Ap/bt/,  ’ L 


(impulsive  realm  beam  bending) 


A = 1.0  for  simply-supported  beams 


•A  = 2.0  for  clamped  beams 


To  e.xpcrinientally  demonstrate  the  validity  of  this  solution,  we  have  plotted  experimental 
data  taken  by  Florence  and  FirthI  ’ ^ I and  comixired  these  data  to  Fq.  (46).  Because  Florence 
and  Firth  used  beams  with  rectangular  cross-sections,  3/,  = Oyhh^/4.  Substituting  for3f, . hit  for 
.1.  and  2t  for  L (they  used  half  spans),  vields: 

7: — = V TT  il  T" ) (impulsive  realm  rectangular  beam  bending)  (47) 

SpOylr  \ V / \ V / 

All  of  the  beams  tested  by  Florence  and  Firth  had  an  C//r  ratio  of  36.  so  this  comparison  is  made 
by  plotting  ij{h\/Npay  ) versus  u’q/C.  All  beams  are  impulsively  loaded  using  sheet  explosive. 

Both  clamped  and  pinned  beams  made  of  2024-T4  aluminum,  606I-T6  aluminum.  1018  cold 
rolled  steel,  and  1018  annealed  steel  are  included  in  this  comparison.  Figure  9 demonstrates  the 
validity  of  Eq.  (47)  and  this  analysis  procedure. 
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Because  we  have  no  data  for  clamped  beams  or  simply-supported  beams  in  the  quasi 
static  realm,  we  will  not  develop  the  equations,  but  wilt  give  the  resulu.  Because  the  clamped 
beam  in  bending  is  8/3  times  Eq.  (42)  for  simply-supported  bending,  we  can  write: 

! j 

1 I2.V*  (quasi-static  realm  beam  bending)  (48) 

I My 

where 


•V  = 1 simply-supported  beam 

V = 2 clamped  beam 


V.  CIRCULAR  PLATE  (BENDING  AND  EXTENSIONAL  BEHAVIOR) 

The  next  solution  that  we  will  evaluate  is  lor  the  residual  mid-point  deformation  in 
uniformly  impulsed  clamped  circular  plates.  This  problem  adds  an  additional  term  to  the  strain 
energv’  expression,  as  both  bending  and  extensional  action  will  be  present.  A possible  deformed 
shape  for  a clamped  plastically  deformed  circular  plate  is: 

3 

/ itr\  * 

w = -^h4-cos~j  (49) 

The  deformed  she;  • is  being  describc-d  by  a radial  (radius  r and  angle  0)  coordinate  system 
with  its  origin  at  the  center  of  the  plate.  Because  of  symmetrx-.  the  deformed  shape  is  inde- 
pendent of  the  angle  0.  Equation  (49)  meets  the  appropriate  boundar>  conditions  for  a clamped 
circular  plate  in  that  at  (he  center  with  r = 0.  »e  is  a ma.ximum  deformation  of  Uo-  and  (du  /dr)  = 

0 while  at  the  edge  of  the  plate  with  r = /?,  u-  = 0 and  (d*e/dr)  = 0-  An  inllection  point  occurs  at 
r - RI2  when  the  curvature  changes  from  negative  to  positive  for  increasing  values  of  r. 

Because  no  change  in  length  occuts  circumferentially,  there  is  no  circumterential  strain  and 
no  circumferential  strain  energy  in  a clamped  circular  plate.  I'he  radial  strain  energy  per  unit 
volume  equals  the  stress  a times  the  strain  c in  any  structural  element.  If  we  assume  that  the 
plate  is  yielding,  the  stre-ss  a must  equal  the  yield  stress  o,  in  a rigid  plastic  constitutive  relation- 
.ship,  and  if  plane  sections  remain  plane,  the  bending  strain  for  small  deformations  is  given  by 
— r(d*  u /dr^ ) where  r is  the  plate  coordinate  perpendicular  to  the  r and  0 coordinates.  1 o com- 
pute the  bending  strain  energy  contribution  L'f,.  create  a circular  differential  torus  of  cireum- 
ferenee  2nr.  thickness  dr.  and  width  dr.  Tire  differential  bending  strain  energy  d( /,  equals 
this  volume  times  o,.  times  r(d'  u/dr- ).  The  total  bending  strain  energy  (■>,  is  obtained  by 
integrating  the  differential  strain  energy  over  the  plate  thickness  and  plate  radius.  Equation  ( 50) 
describes  this  mathematical  procedure. 


/ / {2itrdrAz)(ay){-:-^)  + f f 
L /-*  ,»  \ or  / i./-t  t>  f 

(2ff/-  df  lb)  ( -Oy  ) ^-r  — - J 


Uh 


-h/2  0 


*k/2  H 

s s 

ft/2  «/2 


<50) 


Equution  (SO)  must  be  integrated  in  parts  because  the  sign  of  the  bending  stress  changes 
at  the  inllection  point  which  is  located  at  R\2.  Differentiating  tq.  (49).  substituting  it  into 
tq.  (50).  and  integrating  over  dr  plus  dr  then  yields: 


f'ft  = — 


(51) 


Next  we  estimate  the  extensional  strain  energy  through  a similar  procedure.  I-tlen- 
sional  strains  in  a circular  plate  are  generally  given  by  c = yj\.0  + (du/dr)^  — 1 .0.  Using  the 
binomial  expansion  and  retaining  only  the  first  two  terms  leads  to  the  approximate  expression 
for  extensional  strains,  e ^ 0.5  (du/dr)-.  We  will  use  the  same  ditferential  torus  to  compute 
extensional  strain  energy  as  was  used  to  compute  bending  strain  energy.  The  total  extensional 
strain  energy  l\.  is  obtained  by  integrating  the  differential  volume  times  o»  times  0.5  (dn  /dr  - 
over  plate  thickness  and  plate  radius.  Hquation  (52)  mathematically  describes  this  procedure 
for  obtaining  extensional  strain  energy. 


‘ft/2  R r 1 /du\n 

Differentiating  hq.  (49),  subsiiti'ling  it  into  Hq.  (52),  and  integrating  yields: 


(52) 


L'c  = 


16 


(53) 


The  total  strain  energy  (' is  the  sum  of  (l-'q.(53)l  and  I Eq.  (5 1 )) . or: 


7T*  , 7T’ 

f = — 0,  h-  u „ + — Cyhw  o 
4 16 


(54) 


Tile  kinetic  energy  KE  for  a uniformly  applied  impulsive  load  imparted  to  the  plate  is 
obtained  just  as  it  was  fora  beam,  by  summing  up  the  impulse  squared  divided  by  two  times 
the  incremental  mass  over  the  surface  of  the  entire  plate.  This  summation  leads  to  the  follow- 
ing integration. 


R 

Et:  = / 


0 


/^(2?r>  )^(dr)^ 
2ph{  2wr)  ( dr) 


(5.^) 
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SPECIFIC  IMPULSE 
DENSITY 
YIELD  STRESS 

1.0  PINNED  X 

2.0  CLAMPED 


SYM.  BOUNDARY 


MATERIAL 


CLAMPED 

COLD  ROLLED  STEEL 

■ 

CLAMPED 

2024-T4 

• 

PINNED 

ANNEALED  1018 

▼ 

PINNED 

COLD  ROLLED  1018 

♦ 

PINNED 

2024-T4 

X 

PINNED 

6061-T6 

Data  by  Florence  and  Firth,  Journal  of 

Applied  Mech.  Vol 

. 32.  Sept.  1965,  pp  481 -488 

I Kil'RK  BKAM  HKNOING  IN  Till:  IMPIXSIVK  1 vlADINCI  RK  AI  M 


or: 


ni^K^ 

2p/r 


(56) 


hiuating  I’  (54)|  ti»  AA'  (Hq.  (56)|  yields  the  asymptote  lor  the  impulsive  loading 
r-alm. 


“:\w  >^\h) 


(impulsive  realm  elamped  eircular  plate) 


(57) 


A comparison  between  l:q.  ( 5“)  lor  a unit'orm  impulse  on  a clamped  circular  plate  ami  esperi- 
mental  test  data  can  be  made  usini*  lest  results  by  1‘lorcnce.  I Residual  ivnnanent  mid- 
span delorniations  were  measured  on  both  clamped  circular  bObl-l'b  aluminum  plates  and 
1018-cold  rolled  steel  plates  that  had  been  loaded  unilormly  with  various  layers  of  sheet 
explosive.  The  22  aluminum  data  points.  20  steel  data  points,  and  l;q.  (57)  arc  all  shown  in 
I I}! lire  10.  Once  ag.iin  the  validity  of  this  solution  and  analysis  procedure  arc  subslanliatcd. 

If  there  exists  a svsiematic  error,  it  is  a tendency  lor  the  analytical  curve  to  slijihllv  iinder- 
o-'ihtate  del'ormali<'ns  whenever  w„,'/i  is  lari*e.  This  error  is  probably  caused  by  the  a.ssumcd 
deformed  shape  not  yieldiii)’  the  tnininuun  strain  cncryy.  l or  small  values  of  u ,,///.  the 
anaK  lical  curve  overestimates  deformations.  I'his  err(u  is  created  because  we  assume  that 
deformations  e.xtend  over  the  entire  span  of  the  plate.  When  loads  are  small  ami  defonnations 
small,  the  deformeil  shape  caners  only  a portion  of  the  entire  plate. 


VI.  RECTANGULAR  PLATE  (ADDITION  OF  SHEARING  STRAINS) 

The  final  scriution  that  we  will  develop  is  for  the  plastic  response  of  rectanjrular  plates 
frem  uniformly  applied  transverse  impulses.  This  solution  introduces  a new  complication 
associated  with  a laet  of  radial  symmetry  and  (he  hendiiii!  plus  e.xtensional  shearini;  forces 
that  are  created  thereby.  Fora  elamped-elamped  uniformK  loaded  rectangular  plate,  we  will 
assume  a deformed  shape  yiven  by: 

"'ll  / !i.v\  / m \ 

u---  (\  -Fcos-  ](l  4coSy  j (58) 

where 

.V  ai.d  )'  arc  half  spans 

V and  y are  the  reetanpular  coordinate  system  with  its  oripin  in  the  center  of  the  plate. 

This  a.ssumcd  deformed  shape  meets  the  appropriate  detlection  and  slope  criteria  in  the  middle 
of  the  plate  and  alonp  all  boundaries.  Lines  of  innection  occur  at  y = >72  and  v = ,V.'2.  Ilie 
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strain  eneigy  per  unit  volume  in  a structural  element  under  a biaxial  state  of  stress  is: 
V 


Vol. 


~ J dCjtjf  + 2ojf|.  dcjtj-  + Oyy  1 


(59) 


siruint 


Because  we  have  yielding,  we  will  assume  that  o^k  ~ Oy  and  Oyy  = o,..  but  for  the  shearing 
stress  we  will  use  a Hubcr-Mises-Hencky  distortion  energy  yield  criteria  of  Oxy  - OykJ^- 
The  normal  bending  strains  are  as  in  the  circular  plate  pioblem  with  Cxx  ~ ~c(3*u73v^ ) 
ande,,  = ~r(d*u>/3v* ).  The  bending  shearing  strain  e^,  = 2r(3^iv/3.vdv).  The  bending 
strain  eneigy  Ut,  is  then  given  by: 


HU  X ;■ 
Uh-i  f Azf  d.v  / dr  ov 
0 0 0 


A'  )' 


/ d.v  / di 
0 0 


3^u’\ 

d^J 


dr 


(60) 


- 

y 

u 


i 

r 

r 


The  bending  strain  energy  c'^ntributions  from  the  normal  stress  and  strains  must  be  obtained 
by  four  part  integration  (from  0 to  A'/2  and  0 to  T/2,  from  A'/2  to  A and  0 to  K/2.  from  0 
to  X/2  and  y /2  to  Y.  plus  from  A'/2  to  X and  from  >72  to  K)  because  the  lines  of  inJlection 
cliange  the  sign  of  the  stress.  Differentinting  liq.  (58).  substituting  it  into  Bq.  (60).  and  per- 
forming the  triple  integration  then  yields: 


2 


U'o 


+ -Z,  Oy/l^Wo 


(61) 


The  first  temi  on  the  right  hand  side  of  Bq.  (6i ) is  the  contribution  from  the  normal  bending 
strains,  and  tiie  second  term  is  the  contribution  from  the  bendi;.g  shear  strains. 


Equation  (59)  applies  to  the  extensional  response  as  well  as  bending  response.  The  nor- 
mal extensional  strains  are  as  in  the  circular  plate  problem  with  = 0.5(3u’/3.v)^  and  c,.i-  - 
0.5(3»v/3i  )^ . The  extensional  shearing  strain  = (3>e/3  v)(3u7di  ) as  a first  approximation. 
These  additional  observations  mean  that  the  extensional  strain  energy  (/<•  is  given  by: 


(62) 
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Differentiating  Eq.  (S8),  sutetituting  it  into  Eq.  (62),  and  performing  the  triple  integration 
yields; 


(63) 


The  total  strain  energy  U is  then  the  sum  of  Ut  (Eq.  (63))  and  Ub  (Eq-  (61 )) , or; 


+ Oyhwl 


Y 

X 


4- 


Oyhwl 


(clamped  plate) 


(64) 


The  strain  energy  in  a simply-supported  plate  can  also  be  estimated  using  the  same  pro- 
cedure. with  an  assumed  deformed  shape  described  by; 


ir.v  ity 
w - u 0 cos  — cos  — 


(65) 


The  previous  procedure  gives  as  a value  for  the  strain  energy  (/, 


U ~ Oyh'^  u’o 


Y X A , rr^ 

Ic  '^  Y >/3 


Y X 
X Y 


+ (s.s.  plate) 


(66) 


Notice  that  Eq.  (66)  for  strain  energy  in  a simply-suppoited  plate  is  similar  to  Eq.  (64) 
for  strain  energy  in  a clamped  plate  Tt>e  ur’y  difference  in  these  equations  is  in  the  numeri- 
cal values  of  the  coefficients  which  accompany  each  term  on  the  right  hand  side  of  these 
e.xpressions.  This  observation  means  that  if  we  insert  a coefficient  N that  is  equal  to  1.0  for 
simply-supported  plate  and  equal  to  2.0  for  clamped  plates,  a general  strain  energy  equation 
can  be  written  as  in  Eq.  (67). 


(67) 


The  kinetic  energy  KE  imparted  to  a plate  is  not  dependent  upon  the  deformed  shape. 
For  a uniformly  applied  impulsive  load  imparted  to  a plate,  the  kinetic  energy  is  obtained 
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(as  it  was  For  a circular  plate)  by  summing  up  the  impulse  squared  divided  by  two  tim^  the 
incremental  mass  over  the  suiface  of  the  plate.  The  appropriate  integration  is  given  by: 


J ' 2ph(dx)(dy) 


yv 

(69) 

ph 

Equating  Eq.  (67)  and  Eq.  (69)  finally  yields  a general  rectangular  plate  equ  tion  for 
deformations  caused  by  uniform  impulsive  loads.  The  following  equation  is  this  relationship. 


r 1 

: ^ :r^v-  1 

1 

X 

, 3^-*7r2 

■ 2N 

\fj 

Y 

\h)  16N^ 

Wc  will  demonstrate  the  validity  of  this  solution  by  using  damped  (A'  = 2)  rectangular 
plate  data  reported  by  N.  Jones,  T.  O.  Uran.  and  S.  A.  Tekin.  1*^1  Rectangular  plates  with 
an  aspect  ratio  ( Y/X)  equal  to  1.695  were  loaded  with  sheet  explosive  in  these  experiments. 
Both  hot-rolled  mild  steel  plates  and  606I-T6  aluminum  plates  were  tested  and  can  be  seen 
plotted  in  Figure  I I.  Equation  (69)  is  also  shown  in  Figure  1 1 ; however,  left  and  right  sides 
of  Eq.  (69)  were  multiplied  by  (Y/X)^  to  cast  it  into  the  format  of  their  data.  Excellent 
correlation  appears  in  Figure  1 1 when  Eq.  (70)  is  compared  to  experimental  test  results. 

Unfortunately,  no  experimental  quasi-static  loading  realm  data  exist  on  dynamically 
loaded,  plastically  deformed  plates,  simply-supported  beams,  or  clamped  beams.  The  analyti- 
cal solution  is  easily  developed  for  a plate.  The  work  If  performed  on  either  simply-supported 
or  clamped  plates  is  obtained  by  integrating /ru’  dv  dr  and  equals: 

—2  (A'  - 2> 

h' = - (71  > 

When  this  work  is  equated  to  the  general  strain  energy  expression.  Eq.  (67),  we  obtain  the 
solution  for  the  quasi-static  loading  realm  of  rectangular  plates. 


.J*  r /JCV 


3(A'  - I)4(A^-  A', 


U/J  V'/  2>  [yJV/,  ) 


I^N-  1)41; 

.V-3) 

- 2) 
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VII.  DISCUSSION 


These  procedures  and  the  resulting  formulae  comprise  the  structural  analysis  approach 
that  was  used  by  SwRI  to  control  deformatirms  in  the  design  of  a plastically-yielding,  blast- 
suppressing structure.  This  analysis  procedure  yields  explicit  expressions  for  nondimensional 
deformations  in  terms  of  nondimensional  loading  parameters.  If  other  deformed  shapes  ‘wea* 
to  be  assumed,  the  same  nondimensional  parameters  would  occur;  however,  the  numerical 
coefficients  would  have  other  values.  These  analyses  can  also  be  applied  to  shell  type  struc- 
tures with  no  more  difficulty  than  was  demonstrated  in  the  development  of  the  plate  solu- 
tions. The  approach  is  attractive  because  no  involved  solutions  to  complex  differential  equa- 
tions are  required -only  integrations  of  work,  kinetic  energy,  and  strain  energy  are  needed 
for  a £3ven  structural  element. 

The  illustrations  shown  in  this  paper  have  included  a single-degree-of-freedom  plastic 
model,  a simple  clastic  model,  elastic  cantilever  beams,  plastic  simply-supported  beams, 
plastic  clamped  beams,  plastic  clamped  circular  plates,  plastic  simply-supported  rectangular 
plates,  and  plastic  clamped  rectangular  plates.  Experimental  test  data  obtained  by  other  in- 
vestigators have  been  used  to  demonstrate  the  validity  of  these  solutions,  especially  in  the 
impulsive  loading  lealm.  It  is  interesting  to  note  that  all  plastic  beam  solutions  have  the 
functional  formats 


and 


f * L / w'o\ 

r - y ( — ) (impulsive  realm  any  beam) 

pCy/r  \L  / 


P'-  - . . . . 

r = constant  (quasi-static  realm  any  beam) 

Oyir 


(73) 


(74) 


As  e.vtensional  behavior  is  developed  by  elements  such  as  plates,  these  general  runetional 
formats  arc; 


poyh*  \-Y  h 
and 

o Ij2  '\X'  h / (quasi-static  realm  plates) 

One  could  use  experimental  test  data  that  had  been  nondimensionalized  to  graphically 
present  solutions  following  these  formats. 


(impulsive  realm  plates) 


(75) 
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